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We study the quantum phase diagram and excitation spectrum of the frustrated J1-J2 spin-1/2 
Heisenberg Hamiltonian. A hierarchical mean-field approach, at the heart of which lies the idea 
of identifying relevant degrees of freedom, is developed. Thus, by performing educated, manifestly 
symmetry preserving mean-field approximations, we unveil fundamental properties of the system. 
We then compare various coverings of the square lattice with plaquettes, dimers and other degrees of 
freedom, and show that only the symmetric plaquette covering, which reproduces the original Bravais 
lattice, leads to the known phase diagram. The intermediate quantum paramagnetic phase is shown 
to be a (singlet) plaquette crystal, connected with the neighbouring Neel phase by a continuous 
phase transition. We also introduce fluctuations around the hierarchical mean-fleld solutions, and 
demonstrate that in the paramagnetic phase the ground and flrst excited states are separated by a 
finite gap, which closes in the Neel and columnar phases. Our results suggest that the quantum phase 
transition between Neel and paramagnetic phases can be properly described within the Ginzburg- 
Landau- Wilson paradigm. 
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I. INTRODUCTION 

One of the primary goals of modern condensed mat- 
ter physics is the characterization of strongly correlated 
quantum systems. A large class of such materials is repre- 
sented by frustrated antiferromagnets, which are believed 
to exhibit a variety of novel states of matter at sufficiently 
strong coupling. Growing experimental evidence indi- 
cates that layered materials such as Li2V0(Si, Ge)0^, 
VOMoO^^ and BaCdVO(P04)2^ can be adequately de- 
scribed by an antiferromagnetic Heisenberg model with 
frustrating next- and next-next-nearest neighbor interac- 
tions. As a result, the study of low-dimensional magnets 
and their frustration-driven quantum phase transitions 
have attracted a lot of theoretical attention in the last 
decade^i^. 

A paradigmatic system, illustrating the effects of frus- 
trating couphngs, is the spin-1/2 Heisenberg model on 
a square lattice with competing nearest (Ji), and next- 
nearest (J2) neighbor antiferromagnetic (AF) interac- 
tions (J1-J2 model). Despite numerous analytical and 
numerical efforts, its phase diagram, which exhibits a two 
sublattice Neel AF, quantum paramagnetic, and a four 
sublattice columnar AF states, continues to stir certain 
controversy (for a review of recent achievements, see Ref. 
@). While existence of the Neel-ordered phase at small 
frustration ratio J2/J1, and of the columnar AF state at 
large J2/J1 is widely established, properties of the inter- 
mediate non-magnetic phase, which occurs around the 
maximum frustration value J2/J1 = 1/2, are still under 
debate. Particularly, the correlated nature of the inter- 
mediate state and the kind of quantum phase transition 
separating it from the Neel state, attract most attention. 
Various methods have been recently applied to charac- 
terize the quantum paramagnetic phase, such as Green's 
function Monte Carlo^i^, coupled cluster methods^, se- 
ries expansionsiS and field-theoretical method a^^i^^'^"^ . 



As a result, several possible candidate ground states 
were proposed, namely: spin liquid^, preserving trans- 
lational and rotational symmetries of the lattice, as well 
as various lattice symmetry breaking phases, out of which 
the dimer— li^, and the plaquette resonating valence bond 
phases^ are worth mentioning. 

Not surprisingly, the nature of the quantum phase 
transition separating the Neel-ordered and quantum 
paramagnetic phases is also under scrutiny. The most 
dramatic, and at the same time original, scenario^ is be- 
lieved to violate the Ginzburg-Landau- Wilson paradigm 
of phase transitionsi^ which revolves around the concept 
of an order parameter. Such point of view is based on 
the observation that there are different spontaneously 
broken symmetries in the Neel and quantum paramag- 
netic phases, which thus cannot be connected by a group- 
subgroup relation. The former, of course, breaks the 
SU{2) invariance of the Hamiltonian and lattice trans- 
lational symmetry T— li^, but preserves the four-fold ro- 
tational symmetry of the square C4. On the other hand, 
the paramagnetic phase is known to restore the spin- 
rotational symmetry and is believed to break T and C4, 
due to spontaneous formation of dimers along the links 
of the latticei^ii^. It follows then, that these two phases 
can not be joined by the usual Landau second-order criti- 
cal point. This phase transition can either be of the first- 
order— (the latest coupled cluster calculations^, however, 
seem to rule out this possibility), or represent an exam- 
ple of a second order critical point, which cannot be de- 
scribed in terms of a bulk order parameter, but rather in 
terms of emergent fractional excitations (spinous), which 
become deconfined right at the critical pointi^. 

However, evidence regarding the structure of the non- 
magnetic phase is quite controversial. Indeed, the re- 
sults of spin-wave calculations^, large- iV expansionai^, 
and calculations using the density matrix renormaliza- 
tion group combined with Monte-Carlo simulations^ are 
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believed to indicate the emergence of a dimcr order. 
On the other hand, Monte-Carlo^ and coupled cluster 
calculations^, and analytical results^i seem to support 
the presence of C4 symmetry (plaquette-type ordering) 
in the paramagnetic phase. In the absence of a reliable 
numerical or analytical proof of existence of any partic- 
ular order in the non-magnetic region, there is no ap- 
parent reason to believe in the exotic deconfined quan- 
tum criticality scenario. Although there apparently ex- 
ists numerical evidence^S, at the moment of writing the 
authors are unaware of a local Hamiltonian in space di- 
mensions larger than one, rigorously proven to exhibit 
the type of quantum critical point discussed in Ref. [H. 
Interestingly, it was demonstrated in Ref. Ull that a two- 
dimensional (2D) lattice model can possess a first order 
quantum critical point, which exhibits deconfined excita- 
tions. 

All in all, the complexity of methods used to infer prop- 
erties of the paramagnetic phase and the variety of dif- 
ferent conclusions have created a certain degree of con- 
fusion. Our goal in the present paper is to try to clarify 
some of this controversy by proposing a controlled and 
manifestly symmetry preserving method, geared to com- 
puting ground state properties of the J1-J2 model. Our 
approach is based on the recently proposed systematic 
methodology to investigate the behavior of strongly cou- 
pled systems^^, whose main idea consists of identifying 
relevant degrees of freedom and performing an educated 
approximation, called the hierarchical mean- field (HMF). 
to uncover the phase diagram and other properties of the 
system of interest. In a future work these ideas will be 
coupled to a new, variational with respect to the energy, 
renormalization group approach, which thus adapts to 
the concept of relevant degrees of freedom. 

In the present work we construct HMF approximations 
for the J1-J2 model. The crux of our method is the iden- 
tification of a plaquette (spin cluster 2 x 2 or even larger 
4x4 (superplaquette) symmetry-preserving cluster) as 
the relevant elementary degree of freedom, which cap- 
tures necessary quantum correlations to represent essen- 
tial features of the phase diagram. The importance of 
this degree of freedom was realized only recently in the 
present context^, and somewhat earlier in connection 
with SU{4:) spin-orbital2^, and Hubbard^i models. Be- 
sides being variational, our formalism has the attractive 
feature of preserving fundamental lattice point symme- 
tries and the SU{2) symmetry of the Hamiltonian, by 
utilizing the Schwinger boson-type representation and 
Racah algebra technology. Remarkably, such simple 
mean- field calculation already yields all known results, 
concerning the phase diagram of the J1-J2 model, with a 
good accuracy, namely: existence of a Neel-ordered phase 
with antiferromagnetic wavevector (tt, tt) and spin-wave 
type excitations for J2/J1 < 0.42, a non-magnetic in- 
termediate gapped phase, separated by a second order 
quantum phase transition, and a first order transition 
point, which is characterized by the discontinuous disap- 
pearance of the energy gap and connects the paramag- 
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FIG. 1: A schematic phase diagram of the J1-J2 model, sum- 
marizing our results. In each phase we show spontaneously 
broken (framed symbols) and unbroken symmetries (usual 
symbols). The translational invariance is broken along both 
directions in the Neel and paramagnetic phases, and only 
along the j/-direction in the columnar phase. This fact is 
indicated by the subscripts xy and y after T. Conclusions 
regarding the order of the phase transition, separating Neel 
and plaquette crystal phases, as well as symmetries of various 
phases, are based upon extrapolation of our results towards 
the thermodynamic limit. 



netic state with the columnar antiferromagnetic phase at 
(tt, O) and (0, tt) for J2/J1 > 0.66. 

We emphasize that our investigation primarily focuses 
on the symmetry analysis of the various phases. Out of 
many possible coarse graining scenarios, such as cover- 
ing of the 2D lattice with plaquettes, dimers and crosses, 
only the Ci-symmetry preserving plaquette (or superpla- 
quette) covering (which reproduces the original Bravais 
lattice) displays the correct phase diagram. In particu- 
lar, the intermediate paramagnetic phase is shown to be 
a plaquette crystal, which preserves spin and lattice ro- 
tational symmetries. For all other scenarios, including 
dimerized (bond-ordered) phases, we were unable to re- 
produce all known quantum phase transition points of 
the model. 

We notice that the HMF coarse graining procedure 
leads to an explicit breaking of a particular translational 
symmetry. As a result, one can not draw rigorous conclu- 
sions on the order of the phase transitions, based solely 
on a fixed coarse graining. Nevertheless, it is still possi- 
ble to make some predictions, using a finite-size scaling 
of the relevant degree of freedom towards the thermo- 
dynamic limit, where the effects associated with coarse 
graining should disappear. 

Next two sections are devoted to the formulation of 
the HMF approach. Then, we present results of our cal- 
culations and close the paper with a discussion. Our 
main conclusions are summarized in Fig. [l] which em- 
phasizes symmetry relations between different phases of 
the model. 



II. THE PLAQUETTE DEGREE OF FREEDOM 

We consider the spin-1/2 antiferromagnetic Heisenberg 
model with frustrated next-nearest neighbor interactions 
J2, defined on a 2D bipartite lattice with N sites: 



H 



J: 



(1) 



As mentioned already in the Introduction, we choose the 
plaquette, Fig. [21 as our elementary degree of freedom. 
Then, assuming that N is chosen appropriately, the en- 
tire lattice can be covered with such plaquettes in a sub- 
exponentially^^ (~ 2^) large number of ways. 

Aiming at illustrating the main idea of the method, 
in this section we consider in detail only the symmetric 
covering of the lattice with 2x2 plaquettes, which pre- 
serves the C4 lattice symmetry, see Fig. [Sj although later 
the displaced covering (Fig. [3]), which breaks C4 down 
to C2 (two- fold symmetry axis), and the case of larger 
plaquettes (superplaquettes. Fig. [T51) will be analyzed as 
well. 

It is convenient to take basis the states 



\a) = IhhLM), 



(2) 



where h = Si + S4 and 12 = 82 + S3 are total spins of 
the plaquette diagonals, while L = li + I2 is the total 
spin of the entire plaquette and A/ is its z-component. 
In this basis the Hamiltonian of a single plaquette, 

Ha = Ji {Si + Si) {S2 + S3) + J2 {Si ■Si + S2- Ss) (3) 

is diagonal with eigenvalues 

^hhL [L{L + l)-li {h + 1) - l2{l2 + 1)] + 

+:|[;i(;i + i) + ;2(;2 + i)-3]. (4) 

We note that the basis of Eq. ([2]) is a natural one and 
allows us to explicitly label states with corresponding 
representations of SU{2). 

The next step is to establish how a plaquette couples 
to the rest of the system. In Fig. [3] we show the sym- 
metric plaquette covering of the 2D lattice. In the figure 
the vertices of every non-central plaquette are similarly 
labeled by the numbers 5, 6, 7, 8, and total spins of diago- 
nals are Z3 = S5 + Ss and I4 = Se + Sj. In the uncoupled 
basis matrix elements of the inter-plaquette interaction 




FIG. 2: A single 2x2 plaquette has each vertex occupied 
by a 5 = 1/2 spin. The diagonal spins interact through a 
Heisenberg term of strength J2, while nearest neighbor spins 
interact with strength Ji. 





FIG. 3: Symmetric covering of the 2D lattice with 2x2 pla- 
quettes. Each plaquette is connected to 4 nearest and 4 next- 
nearest neighbors. 



(5) 



LM 



X {L[M[L'2M^\L[L'2LM){LiMiL2M2\LiL2LM), 

where cr = 1 (c = 2) corresponds to the nearest (next- 
nearest) neighbor interaction, Li,L2 {L'i,L'2) represent 
initial (final) angular momenta of the two plaquettes and 



L = Li + L2 is their total angular momentum. In this 
equation we have introduced the notations Ai = {lil2Li\, 
A2 = {'3^4^21 and ai = {\iMi\, and similarly for the 
primed indices. Because each plaquette has 4 nearest 
neighbors and 4 next-nearest neighbors (see Fig. [3]), the 
symmetrized next-nearest neighbor interaction may be 
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written as: 

(A'iAi,LM|iJ?JAiA2,LM) ^ J2pf;g(i)x (6) 

while the symmetrized nearest neighbor plaquettc inter- 
action has the form: 

{X[X'^,LM\HU>^,X2,LM) - Jipf;g(L)x (7) 

+ 2{X[X'2,LM\Hl^\XiX2,LM), 

In Eqs. del) and the symbols S"^'^ = (A'H^nljA) denote 
reduced matrix elements of the n-th spin operator, and: 

n^'i^'2(r \ — ^ ( -WL+Li.+Li \ A L'2 L 1 
Pl.lS^) - 4(-l) \ L2 ii 1 J ' 

where {• • • } are Wigner &j symbols (or Racah 
coefficients)^^. 

Let us now identify the plaquette degree of freedom 
with a Schwinger boson which creates a specific state of 
the plaquette. Then, the Hamiltonian of Eq. ^ in the 
plaquette basis can be expressed as: 

H=^^allal^a + (^hit ) I^a^ 7,1^ 7»ai 7ja2 + 
i,a (ij) 

+ (^i^it) llZ ^]<^: ^^''^ ' 

where the operator creates a boson on site i of the 
plaquette lattice (which contains A^q ~ N /A sites) in the 
state, denoted by an index a, running through the entire 
single-plaquette Hilbert space (of dimension 2^ = 16) 
and the summation is performed over doubly repeated 
dummy indices. The unphysical states are eliminated by 
enforcing the local constraint '^allalia — 1- what 
follows, we impose periodic boundary conditions on the 
plaquette lattice. 

The bosonic operators define the hierarchical 
language^^ for our problem. It will be used in the next 
section, where we develop an approximation scheme for 
diagonalizing the Hamiltonian of Eq. ([5]) . 

III. HIERARCHICAL MEAN-FIELD 
APPROXIMATION 

As it follows from Eq. ([4]), the lowest single-plaquette 
state has the energy eiioo/4 = — Ji/2 -I- J2/8 per spin, 
which, when J2 = 0, gives only the energy of a classical 
2D antiferromagnet. Thus, it is necessary to take into 
account the interaction term in Eq. ([5]). 

The HMF approximation is a mean-field approach, 
performed on the relevant degrees of freedom. In the 



present section we discuss only the simplest one - a 
Hatree-Fock like (HE) approximation. A possible way 
to include fluctuation corrections is presented in Ap- 
pendix B. The HE approximation introduces the mix- 
ing of single-plaquette states which minimizes the total 
energy of the system and is based on a canonical trans- 
formation among the bosons, which we will restrict to be 
uniform (plaquette independent): 

l^a=Ra'^^n■ (9) 

The real matrix R satisfies canonical orthogonality and 
completeness relations: 

nn nn' r . nri nn r 

A translationally invariant variational ansatz for the 
ground state (vacuum) is a boson condensate in the low- 
est HE single-particle energy state (n = 0): 

|HF)=nrlolO>, (10) 

i 

and since it has one boson per plaquette, there is no 
need to impose the Schwinger boson constraint in the 
calculation. 

Minimizing the total energy with respect to R, we ar- 
rive at the self-consistent equation: 

LSaa' +Y.z4H^'XZ,«]K' = enK, (H) 

where zi = Z2 = 4 are the nearest- and next-nearest 
coordination numbers. The ground state energy (GSE) 
per spin is then given by the expression: 

^ a ^ 

with eo being the lowest eigenvalue of Eq. (fTT|) . 

Another fundamental quantity to compute is the po- 
larization of spins within a plaquette: 

where n = 1,...,4 is the spin index, and the matrix 
elements (determined from the Wigner-Eckart theorem) 
are: 

(^fJa'a = {l'AL'M'\S:,\hhLM) = (13) 

( ^^L+L'+l. (10£M|1££'M) ,„„^„,^ 
= (-1) ^ ^ dMM' ^ ^ {IJ^LWSnWhhL}. 

This enables us to define the staggered and coUinear 
(along X and y axes) magnetizations: 

M,tag - (1/4) imSt + SI - - |HE) ; (14) 
Afcoi(x, y) = (1/4)(HE|5^ - 5| + 5|,3 - Sl^lRF). 

Notice the extreme simplicity of the HME approxima- 
tion. The reason why it is able to realize meaningful 



FIG. 4: The displaced plaquette covering. Notice that the C4 
symmetry is broken down to C2. 



FIG. 5: Connectivity of the dimer lattice for symmetric (right 
panel) and displaced dimer coverings (left panel). The rota- 
tional C4 symmetry is lowered to C2 in both cases. 



results is that the plaquette degree of freedom seems to 
contain the main correlations defining the physics behind 
the Hamiltonian of Eqs. ([T]), and (HJ. To avoid confu- 
sion, we emphasize that the HF approximation and the 
fluctuation theory of Appendix B arc physically (and ob- 
viously mathematically) different from the spin-wave or 
canonical Schwinger-Wigner boson mean-field approach 
to spin systems^. In particular, we make no assumption 
about the underlying ground state, thus allowing for an 
interplay of various quantum phases. Moreover, it will 
be demonstrated, that the collective excitation spectra 
in each phase consistently reflect spontaneously broken 
symmetries, unlike the usual Schwinger boson caso2S, in 
which one obtains gapped excitations. 



FIG. 6: Covering of the lattice with crosses - arrays of five 
spins. Since one cross cannot form a singlet, it is necessary to 
double the unit cell, as indicated by the gray shading. This 
choice of a degree of freedom clearly preserves the C4 symme- 
try, but the resulting lattice breaks it. 



IV. GROUND STATE PROPERTIES AND 
EXCITATION SPECTRUM OF THE MODEL 

Our choice of the plaquette as an elementary degree 
of freedom remains unjustified at this point. In order to 
show its relevance we applied the analysis of two previ- 
ous sections to several other coarse grainings (besides the 
symmetric plaquette covering, case (a), shown in Fig. [3]): 

(b) superplaquette (spin cluster 4x4) degree of freedom, 
covering the lattice in such a way that C4 is preserved 
(see Appendix A for details); (c) displaced plaquette cov- 
ering of the lattice. Fig. |4l (d) symmetric and displaced 
dimer coverings, shown respectively in right and left pan- 
els of Fig. [5l (e) cross degree of freedom, Fig. [6l One 
should observe that symmetries of the original Bravais 
lattice are preserved only in cases (a) and (b). In cases 

(c) and (d) the lattice rotational symmetry C4 is lowered 
to C2. Case (e) is special in the sense that an isolated 
degree of freedom does not possess a singlet ground state. 
The information about a particular configuration is en- 
coded in matrix elements of -ffiJ^tj whose calculation is 
elementary. Other equations, presented in Sees. II and 
III, retain their form. 

For each of the above cases we iteratively solve Eq. (fTTj) 
and compute the GSE (|12|) . and staggered and coUinear 
magnetizations (|14p . The main message, which wc would 
like to convey in this section is that only the plaquette 
degree of freedom (of any size) is relevant for constructing 
the phase diagram of the Hamiltonian of Eq. ([T|) . 



A. Symmetry preserving plaquette configurations 

Let us focus first on cases (a) and (b), i.e. symmetry- 
preserving coverings of the lattice with plaquette. Fig. 
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FIG. 7; Ground state energy per spin computed at the HF 
level for the 2x2 and 4x4 plaquette elementary degrees of 
freedom. The inset shows finite-size scaling in the AF phase 
at J2 = 0. 



[3l and superplaquette degrees of freedom. The resulting 
GSE as a function of J2/J1 is shown in Fig. [71 One im- 
mediately observes a level-crossing at w 0.67Ji, indi- 
cating the first-order transition and a second-order quan- 
tum critical point at JJ^ ~ 0.42Ji, which is supported 
by a jump of the second order derivative (fEo/dJ2, 
Fig. [51 Both Neel and columnar phases are charac- 
terized by spontaneously broken SU{2) symmetry. The 
former exhibits a nonvanishing staggered magnetization, 
A/stag, while the latter has nonzero coUinear magnetiza- 
tion along the x-direction, Mcoi(a;). Both order parame- 
ters become zero in the paramagnetic phase, suggesting 
that SU{2) is restored. These results are summarized in 
Fig. [51 from which it also follows that the phase tran- 
sition at is continuous, while corresponds to a 
first-order transition point. We remind, in this connec- 
tion, that our approach does not explicitly break the spin 
rotational symmetry, thus allowing for the treatment of 
competing ground states. 

As expected, considering a larger elementary degree of 
freedom - superplaquette - leads to a significant improve- 
ment of the GSE and reduction of the magnetization, due 
to larger quantum fluctuations. The finite-size scaling 
(insets in Figs. [7land[9|), using these two sizes (2x2 
and 4x4), indicates that i?o(J2 = 0)/7V -0.64Ji and 
Mz{J2 = 0) ^ 0.39 in the thermodynamic limit, a sat- 
isfying result for a HF approximation, which completely 
ignores fluctuations (these numbers should be compared 
with well-known results of Monte-Carlo simulations^!: 
Eq/N k -0.67 Ji and w 0.31). 

Next, we discuss in more detail symmetry properties of 
the various phases in Figs. [3 [3 At all values of J2 < J|^, 
the lattice translational symmetry T is brokeni^, but the 
rotational C4 symmetry is preserved. For J2 < this 
corresponds to a Ncel-type long-range order with sponta- 
neously broken SU{2). At large values J2 > JJ^, we ob- 



FIG. 8: Second-order derivative SEo/dJl for the 2x2 and 4x 
4 plaquette degrees of freedom. The discontinuity at J2/ Ji ~ 
0.42 is indicative of a second-order quantum phase transition. 
In the inset we present finite-size scaling for the jump g = 

(Jid2i5o/iVdj|)J+°. 



0.5 




J2/J1 

FIG. 9: Staggered magnetization, Matag, for J2 < J^^ and 
coUinear magnetization along the i-direction, Mcoi(2;), for 
J2 > J2^ (for the 2x2 and 4x4 plaquette degrees of free- 
dom) , computed at the HF level. Notice the continuous phase 
transition at J2/J1 ~ 0.42 and a first order transition at 
J2/J1 ^ 0.68 (2 X 2) and J2/J1 « 0.66 (4x4). The inset 
shows finite-size scaling of Mz at J2 = 0. 



serve the columnar ordering, which spontaneously breaks 
C4 down to C2, and SU{2), but partially (i.e., along 
one direction) restores the lattice translational symme- 
try. We present a more detailed discussion of the spatial 
symmetries later in this section. In the intermediate re- 
gion J2 G (^2^, J2^) the spin 5/7(2) rotational symmetry 
is restored. In this paramagnetic phase the ground state 
wavefunction is a tensor product of individual plaque- 
tte ground states (with quantum numbers li = I2 = 1, 
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L = M = 0): 
|1100) 



2V3 



2(|TlT4i2i3) + Uli4T2T3»- 
-(|Tli4T2i3)+|Tli4i2T3>+ (15) 
+ UlT4T2i3) +1 ilT4i2T3)) . 



This ground state necessarily breaks the lattice transla- 
tional symmetry, but preserves C4. In fact, the param- 
agnetic region on the phase diagram of Figs. [7] and [9] is 
a trivial plaquette crystal: a set of non-interacting pla- 
quettes, because the expectation value of the plaquette 
interaction (see Eq. ([5])) in the singlet state fl- llnool^) 
vanishes. An analogous situation is realized when the 
superplaquette is chosen as an elementary degree of free- 
dom: the paramagnetic phase is a crystal of superplaque- 
ttes. It is interesting to note that in Ref. a "plaquette 
resonating- valence-bond state" , exactly equal to ([15]) , has 
been proposed. However, later— the intermediate phase 
was argued to be a spin liquid, i.e., a state that preserves 
the lattice translational symmetry. 

In order to learn about spatial symmetries in various 
phases, we compare magnitudes of the several lattice 
symmetry-breaking observables proposed in the litera- 
ture. We consider the following three, introduced in Ref. 

(in the notation of that paper): 

-^1 =ly y^X^^y Sx,ySx + l,y] 

F2 ^ Sx,y{Sx+l.y - Sx.y+l)] 

x,y 

Fi =^ ^--V [{-^TS., + ,.y + (-l)^S.,,+ l] , 

x,y 

where indices x, y specify a spin in the 2D lattice. The 
operator probes the plaquette ordering, which pre- 
serves the lattice rotational symmetry, while Fi and F2 
correspond to the columnar ordering. We note, however, 
that Fi is already non-zero for an isolated plaquette (or 
superplaquette). These functions can be combined in 
the complex "order parameter", introduced in Ref. [Tsl . 
Here we show details of the calculation of functions ^1,2,4 
for the plaquette degree of freedom, case (a), and only 
present the result for F/^ for the superplaquette case (b). 
In the plaquette representation the above operators are 
written as: 



^ ('S'2;i,j5^5;i-l-l,i + 'S'4;ijS'7;i-|-i j)] ; 
^2.4 =TF Etl'^l^^.J- T Si■,^.,){S2■,^., T 53;.,,)± (16) 



± ('S'2:i,j'S'5;i-|-i J- -I- S'4;i j-5'7;i+i 
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FIG. 10: HF ground state expectation values of the 
symmetry-breaking perturbations, given by Eq. (|16p . plotted 
as functions of J2/J1 for the plaquette and superplaquette 
degrees of freedom. Due to the unbroken C4 symmetry in the 
Neel and paramagnetic phases, values of F4, are twice larger 
than the corresponding values of F\, except in the columnar 
phase, where C4 is spontaneously broken down to C2. The 
inset shows finite-size scaling for Fi for three values of J2 / Ji : 
(circles), 0.504 (triangles) and 0.997 (rhombs). 




FIG. 11: The two lowest excitation energies taken at the cen- 
ter of the plaquette Brillouin zone. The main panel shows 
the self-consistent solution to Bogoliubov's equations, while 
the inset corresponds to the time-dependent Gross-Pitaevskii 
equation (weak coupling). Since wavefunctions of collective 
excitations in the Neel and columnar phases have different 
symmetries, there are level crossings in the non-magnetic 
phase (cusps in the plot). 



In this equation the indices i, j are coordinates of a pla- 
quette in the lattice. 

Expectation values of the functions Eq. (fT6)l in the HF 
ground state are shown in Fig. 1101 Both phase transi- 
tion points and Jf^ are clearly seen from this plot. 
All functions change continuously across the second-order 
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critical point J!^^ and jump at the first-order transition 
point J|^. Except in the columnar phase the values of 
F4 are everywhere exactly twice larger than those of Fi, 
which is an indication of the unbroken four-fold rotational 
symmetry of the lattice in these regions. In the columnar 
phase, on the other hand, this symmetry is broken and 
the above relation docs not hold. While in the Ncel and 
columnar phases nonlocal terms in Eq. (|16p arc impor- 
tant, in the paramagnetic state the only contribution to 
either expectation value comes from isolated plaquettes 
(local terms in Eq. HH])), or superplaquettes. This obser- 
vation is consistent with properties of the ground state in 
the non-magnetic phase, discussed earlier in this section. 

As mentioned in the Introduction, any choice of degree 
of freedom breaks explicitly the lattice translational sym- 
metry r, with the result that links in the lattice become 
inequivalent. Indeed, the functions of Eq. (jl6p . defined 
on the links, have non-zero values even in the AF phase 
at J2 — 0. However, this effect vanishes in the ther- 
modynamic limit (i.e., as the size of the degree of free- 
dom is increased). The finite-size scaling for F4{J2/Ji) 
is presented in the inset to Fig. [TUl The three extrap- 
olated values: ^4(0) = -0.075, F4 (0.504) = -0.245 and 
^4 (0.997) = -0.004, suggest that the "hnk-wise" transla- 
tional invariance is restored in the thermodynamic limit 
in the Ncel and columnar phases, but not in the plaquc- 
ttc crystal phase. Moreover, the value of the jump g = 

'^pfdf° I j^i-o' extrapolated to the thermodynamic limit 
(see inset to Fig. |8]) remains finite: g{Ji^} = —1.006. In 
other words, these results imply that the critical point 
corresponds to the usual Landau second-order phase 
transition. 
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FIG. 12: HF ground state energy per spin for the displaced 
plaquette and cross (upper panel), and dimer (lower panel) 
coverings. 



B. Excitations in the plaquette crystal phase 

Until now we have considered only ground-state prop- 
erties of the model Eq. (jT]). However, low- lying excited 
states are also of considerable interest. In particular, 
the paramagnetic phase is known to have gapped exci- 
tations, while Ncel and columnar phases exhibit Gold- 
stone modes. Thus, the phase transition points J^^ and 
must be accompanied by the opening of a gap in 
the excitation spectrum: the former in a continuous and 
the latter in a discontinuous fashion. In Appendix B 
we present a particular method to obtain the collective 
spectrum of the system. The main idea of this approx- 
imation is borrowed from the Bogoliubov-Fetter theory 
of superfluidity--'^. Namely, assume that on each plaque- 
tte the majority of Schwinger bosons form a condensate 
in an appropriately chosen lowest energy state and ne- 
glect fluctuations in the number of condensed particles. 
We note, however, that due to the Schwinger boson con- 
straint, this quantity has the meaning of a probability 
to find a given plaquette in the lowest energy HF state, 
rather than the number of particles. Nevertheless, we 
will call it the condensate fraction hq, which, in principle. 



should be determined self-consistently, and is a measure 
of the applicability of the entire approximation: it should 
satisfy the inequality |no — 1| <^ 1. Once the condensa- 
tion part is separated from jia, what remains describes 
fluctuation corrections to the HF ground state. These 
fluctuations have rather strong effects near the phase 
transition points, leading to the modification of to 
a A-point and its considerable shift. The value of 
also changes, but much less significantly. These facts im- 
ply that our approximation breaks down near the phase 
transition points. Indeed, in Appendix B it is shown 
that close to the transition, the condensate is strongly 
suppressed. However, deep in each phase uq ^ 0.9, thus 
allowing us to draw conclusions about general properties 
of the collective spectrum. 

The complete summary of the results is given in Ap- 
pendix B, here we only present the most interesting one: 
The gap in the excitation spectrum as a function of 
J2/J1. Although we focus only on case (a) (2x2 plaque- 
tte), the superplaquette degree of freedom can be con- 
sidered in a similar manner. The gap always occurs in 
momentum space at fc = 0, which reflects translational 
invariance of the plaquette lattice. Below wc focus only 
on this point in the plaquette Brillouin zone. In fact, 
there are 16 — 1 = 15 collective branches and only some 
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of them become gapless in the phases with spontaneously 
broken SU{2). However, in the paramagnetic phase all 
branches develop a gap. In Fig. [Tl] we show the en- 
ergy gap A( J2) = uj{k = 0, J2) for the two lowest ex- 
citation branches in a system of 100 x 100 plaquettes, 
which approximates well the thermodynamic limit. The 
main panel shows results of the self-consistent solution of 
Bogoliubov's equations. The inset compares it with the 
solution of the time-dependent Gross-Pitaeveskii equa- 
tion (pS)) . which corresponds to the weak-coupling ap- 
proximation. In the Neel and columnar phases there are 
two spin wave-type Goldstone modes, both of which ac- 
quire a gap in the paramagnetic phase, at and J|^. 
However, as it follows from Fig. [Til positions of these 
points change from their HF values to: sa 0.33 Ji and 
J|2 w 0.65Ji. The critical point was obtained by 
extrapolation of the staggered magnetization curve (Fig. 
[THll to zero, while the first-order point by extrapolat- 
ing the two GSE curves in Fig. [T4]until intersection. The 
single-plaquette physical picture, discussed previously in 
connection with the paramagnetic phase, remains valid, 
e.g. the condensation occurs again in the plaquette state 
1 1100). Using this observation and symmetries of the 

matrix elements (^Hf^^)"^!^ , one can rigorously show the 
existence of a gap in the non-magnetic region. In fact, we 
can say, that it is a property of our HMF approximation, 
rather than a numerical evidence. 



C. Other degrees of freedom 

Finally, we comment on the results for cases (c)-(e), 
which, contrary to the configurations considered before, 
explicitly break lattice rotational symmetry, see Figs. [4]- 
ini The corresponding GSEs are shown in Fig. [T^l In con- 
trast to the previously considered scenarios, these cases 
give qualitatively wrong phase diagrams. Indeed, if we 
cover the lattice with displaced plaquettes or crosses, 
there exists no classical spin configuration, which gives 
the long-range Neel order. On the other hand, such con- 
figuration exists for the columnar state. For the displaced 
plaquette covering the low- J2 phase J2 < is an SU (2) 
singlet, and spatially is a set of non-interacting plaquettes 
(notice the coincidence of 2 x 2-plaquette energies in the 
paramagnetic phases of Figs. [71 and fT2| . Thus, the phase 
transition to the columnar state is of the first-order. For 
the cross covering, on the other hand, SU{2) is explic- 
itly broken for all values of J2, but since the columnar 
phase partially restores the lattice translational invari- 
ance, it is again separated from the non-magnetic state 
by a first-order phase transition point. The two dimer 
configurations, case (d), are complementary to each other 
in the sense that one of them has only the classical Neel 
state and another - only columnar phase. It follows that 
these configurations can have only one second-order crit- 
ical point at which SU{2) is restored and other symme- 
tries remain broken. As a result, one obtains the phase 
diagram, shown in the lower panel of Fig. 1121 which is 



invariant under reflection in the plane J2 = J|^. These 
observations imply that the coarse graining prescriptions 
(c)-(e) are probably a bad starting point for any approx- 
imation scheme. 



V. DISCUSSION 

In this section we would like to put our main results in 
perspective by making several summarizing remarks. It 
should be emphasized, that the discussion below is based 
on our finite-size scaling results. 

First of all, from our analysis it follows that dimer 
(bond) order is always unfavorable in the non-magnetic 
phase. Notice that even when the plaquette coverings 
were considered such an order did not occur, although 
spontaneous dimerization was not explicitly prohibited. 
Instead, the quantum paramagnetic phase prefers to 
preserve the lattice rotational symmetry, which makes 
the phase transition separating it from the Neel phase 
fit perfectly well within the Ginzburg-Landau- Wilson 
paradigm. The data presented for the staggered magneti- 
zation. Fig. [HI and symmetry-breaking observables. Fig. 
[TOl indicate that the symmetry group of the Neel state is 
a subgroup of the symmetry group in the paramagnetic 
phase, as both phases break T and preserve C4, but the 
latter also preserves SU(2). On the contrary, there is no 
such group-subgroup relation between the paramagnetic 
and columnar antiferromagnetic phases. Consequently, 
the transition between these two states is first order. 
These observations are summarized in Fig. [TJ Indeed, 
starting from the known symmetry in the Neel state and 
assuming validity of the Landau theory, one can unam- 
biguously rule out dimerized structures in the paramag- 
netic phase, since they break lattice rotational symme- 
try. Therefore, our results do not favor the scenario of 
deconfined quantum criticality, advocated in Refs. [^Iisl. 
As already discussed in the previous section, our method 
explicitly breaks a particular lattice translational sym- 
metry: The ground state in Eq. ([10|) at J2 = (AF 
phase), is not invariant under a [11] lattice translation. 
One way to cure this problem is to consider variational 
wavefunctions of the "resonating plaquette" type: 

1*) = (l + rn)|HF), 

which restore that symmetry, with Tn the translation op- 
erator along the [11] direction in the lattice. This state 
describes two resonating plaquette configurations, shifted 
with respect to each other along [11]. Results of calcula- 
tions using this wavefunction for systems up to 6 x 6 spins 
indicate that for J2 < the ground state has long- 
range Neel order and is paramagnetic for J2 € (J|"'^, J^^)- 
The intermediate phase has a plaquette crystal order, 
but with partially restored translational invariance. The 
phase transition at J2^ is still of the second-order, which 
is not surprising, as it can be described solely in terms of 
the SU (2) order parameter. However, based on these sys- 
tem sizes, we can not definitively conclude whether this 
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phase transition remains of the second-order or becomes 
weakly first-order in the thermodynamic Umit. 

Next, wc observe that despite profound differences be- 
tween the 2D and ID equivalent J1-J2 models, their 
non-magnetic phases present some similarities. The 
one-dimensional model is known to be quasi-exactly 
solvable^ at the point J2 = 0.5 Ji and exhibits a param- 
agnetic ground state with short-range correlations for J2 
above the critical valuers J| w 0.24Ji (however, due to 
the peculiar physics in one dimension, the critical point 
J 2 is an essential singularity and, therefore, not obvi- 
ously accessible for the HMF approximation of the type 
presented here). In this non-magnetic region the ground 
state is doubly degenerate, corresponding to two possi- 
ble coverings of a ID lattice with dimers, in accordance 
with the Lieb-Schultz-Mattis theorem^l. Unfortunately, 
a finite-size scaling calculation for the gap between the 
lowest and first excited energy levels, based on exact di- 
agonalization of the 2x2 and 4x4 clusters with peri- 
odic boundary conditions, does not provide a definitive 
answer to the question on whether the ground state of 
the 2D J1-J2 model becomes degenerate in the region 
J 2^ < ^2 < ■ This is indeed what one would ex- 
pect on the basis of a generalization of the Lieb-Schultz- 
Mattis theorem to higher space dimensions (see, e.g. Ref. 
[12). At the HF level, it is true that different plaque- 
tte coverings of the lattice have the same energy (simply 
because each plaquette is in its singlet ground state). 
However, the total number of such configurations grows 
sub- exponentially ^ 2^ , which should be contrasted 
with the dimer covering problem, where this number is 
known^ to be exponentially large. Based on this distinc- 
tion, one may speculate that if our plaquette picture is 
valid, there are not enough different plaquette configura- 
tions for the paramagnetic phase to become a spin liquid 
(i.e., a resonating plaquette state). This statement, cer- 
tainly, requires a separate investigation. 

Finally, we emphasize that the main goal of this work is 
to investigate the fundamental symmetries of the phases 
exhibited by the J1-J2 model. Although the energies 
presented for the 2x2 and 4x4 plaquette cases are 
different from those, obtained by more sophisticated nu- 
merical methods, they can be systematically improved 
by considering correlated trial wavefunctions or by using 
more complex methods, which build upon the results re- 
ported here. However, we expect that symmetry-wise our 
conclusions will remain unchanged. One of such methods 
suited for computing the phase diagram of the model H]), 
which received significantly less attention, amounts to 
applying the Wilson renormalization group proceduro^ 
and the density matrix renormalization group (DMRG) 
method^. Application of the latter faces serious difh- 
culties in 2D (see Ref. HH for a related discussion) and 
requires a mapping of the 2D system onto a chain, which 
introduces a certain bias to the final results^. While for- 
mulation of a practical and efficient DMRG approach to 
2D systems is yet to be developed, we note that our re- 
sults may provide a useful and rational initial input for 
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FIG. 13: The 4x4 superplaquette degree of freedom. Each 
spin carries two indices: a 2 x 2 plaquette number and a 
coordinate within this plaquette. 



such an algorithm. 



VI. CONCLUSION 

In summary, we analyzed the phase diagram of the 2D 
J1-J2 model on a square lattice, focusing on symmetries 
of the various phases. 

We showed that in this model the hierarchical 
language^^ is defined by identifying the plaquette as a 
relevant degree of freedom. Using an unbiased and mani- 
festly symmetry-conserving mean-field approach we com- 
pared it with several other possible candidates: dimer 
and cross degrees of freedom, as well as different ways to 
cover a 2D lattice with plaquettes. Our results indicate 
that the plaquette (and superplaquette) covering, which 
preserves the lattice rotational symmetry has the best en- 
ergy among considered coarse graining scenarios and it is 
the only one to reproduce known facts, such as the inter- 
mediate phase with gapped excitations, concerning the 
phase diagram of the model, while other configurations 
fail to exhibit all quantum phase transition points. 

Consistent with the previous work on the subject, we 
found the quantum paramagnetic phase in the interval 
0.42 < J2/J1 < 0.66. Main controversies revolve around 
the nature of this intermediate non-magnetic phase and 
of the quantum critical point separating it from the Neel 
- ordered state. We found that the paramagnetic phase 
is a plaquette crystal, preserving both lattice and spin 
rotational symmetries. Extrapolation of our numerical 
results to the thermodynamic limit suggests that the 
Ginzburg-Landau- Wilson paradigm of phase transitions 
is perfectly applicable in this case. Indeed, within the 
HMF, there is a group-subgroup relation between sym- 
metries of the non-magnetic and Neel phases, which are 
thus separated by a second-order phase transition. On 
the contrary, such relation does not exist between the pla- 
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quettc crystal and columnar antiferromagnetic phases, so 
the corresponding transition is of the first order. Our pla- 
quette crystal is quite different from the usually proposed 
dimerized (bond-ordered) phases in this non-magnetic re- 
gion. 

We also proposed a way to include fluctuations around 
the HF ground state and showed that properties of the 
collective excitation spectrum are consistent with the 
overall picture of spontaneously broken symmetries. In 
particular, it was demonstrated that the quantum param- 
agnetic state is characterized by a finite gap in the excita- 
tion spectrum, which vanishes in the Neel and columnar 
phases, producing a doubly degenerate Goldstone mode. 

Although currently there exists no known material 
whose ground state realizes the paramagnetic phase of 
the J1-J2 model, in the future momentum-resolved mea- 
surements, such as neutron diffraction methods, can be 
used to identify the plaquette crystal phase of the type 
proposed here. Its experimental signature will be the 
unbroken four-fold lattice rotational symmetry on both 
sides of the second-order phase transition at the critical 
point J^^. 

We acknowledge fruitful discussions with CD. Batista, 
C. Esebbag and E. Fradkin. This work was supported in 
part by the Spanish MEG under grant No. FIS2006- 
12783-C03-01. 



VII. APPENDIX A: SUPERPLAQUETTE 
DEGREE OF FREEDOM 

Here we present details of the HF calculation which 
uses the 4x4 superplaquette, shown in Fig. [T31 as an 



elementary degree of freedom. It turns out that the full 
angular momentum basis is inconvenient, thus we use 
the 2x2 plaquette product states in order to perform 
the mean-field calculations. Each spin is characterized 
by two indices: the plaquette number i = 1, ... ,4 and 
an index n = 1 ... 4, which specifies a vertex in the pla- 
quette. The singlet sector of the superplaquette Hilbert 
space is spanned by the states: 



4 

\A)=]X\a.), 

i=l 



where prime indicates the constraint X]i=i ~ ^- Us- 
ing these states, we can write down matrix elements like 



{a[ . . . a'^\SinSjn' \ai . . . 04) 



in the compact form: 

{jV\S,„S,n'\A) ^ M"^^ l[S,,,, (17) 

with the symmetric matrices (o'nm)°'°^ = (''^"™)a^a' ~ 
{(^nm)l'2' defined as: 



K,M 



% 2 f }(A:||5„||A.)(A;||5„,||A,), 

The Hamiltonian of a single superplaquette consists of two parts: a diagonal one, involving only 2x2 plaquette 
contributions and a non-diagonal part, which accounts for the plaquette interactions. The former is written down 
straightforwardly as a matrix: 

i i 

where is the plaquette self-energy, Eq. The non-diagonal part has the operator form: 

^^nd = Jl[(5'i4 + 5'4i)(523 + 5^32) + S'i2S'21 + Sr^Sj,i + 5^34543 + 5'24'S'42] + 

+ J2 [5^125^23 + 'S'l4'S'21 + 5^235^42 + 'S'4iS'24 + 5^345^41 + S^32'S'43+ (19) 
+ 'S'3iS'i4 + S13S32 + 5'i45'4i + 5^325^23] ■ 

Let us now proceed with interaction terms in the J1-J2 Hamiltonian ([1]). Each superplaquette has 4 nearest and 
4 next-nearest neighbors. Within each neighboring superplaquette we enumerate 2x2 plaquettes by the indices 5, 
6, 7, 8, so that 1— *5, 2-^-6, 3— >7 and 4^8. Enumeration of vertices within a plaquette stays the same. In this 
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manner wc have the symmetrized nearest neighbor 



-ffl — — [Sll (S'73 + S'62) + S'i2S'74 + Si3S(j4 



83 



S31SS2 + + S'33(S'84 + '§'51) + 

+ S22{Ss4 + S51) + 8433^1 + S42SJI + S44{Sj3 + 'S'62)]- 
Sl2{Sj3 + Ss3) + 813(3^2 + Ss2) + S2l{S74 + Ss4) 
^83) + 824(851 + S71) + S'3i(5'84 + 854) + S33{Ss2 + '§'52) + 



'S'22('S'53 
834(851 + Sei) + 842(373 



S53) + S43(Sq2 + S52) 
I 



(20) 



and next-nearest neighbor 
J2 



H2 



['S'llS^84 + S22S73 + S'33S^62 + S'44S'5i] (21) 



superplaquette interactions. Using Eq. (|17p . one can 
easily construct matrix elements of the operators (|19|) - 
(PT|) . which are required to obtain the HF equation of 
the type pT|) . 

Having computed the singie-superplaquette ground 
state wavefunction R\ = R^aia2a3a4,i we can use it to 
determine the spin polarizations: 

(HF|5i„|HF) = (S',';)^,^^^i?°/^a2a3a4-^aia2a304! 
(HF|S'|„|HF) = (S'fi)^^^^i?Qja2a3a^-Raia2a3a4' 

where {S^)^,^ is given by Eq. 



VIII. APPENDIX B: FLUCTUATION 
CORRECTIONS - SUPERFLUID MEAN-FIELD 

In this Appendix we extend the analysis of Sec. Ill 
by considering fluctuations around the HF ground state. 
While not unique, a natural way to achieve this goal is 
to perform a superfluid-type mean-field approximation. 
As a result one can obtain the collective spectrum and 
corrections to the GSE and magnetization. Of primary 
interest is, of course, the energy gap in the excitation 
spectrum. 

The structure of the superfluid mean-field is similar to 
the Fetter-Bogoliubov approach to inhomogeneous Bose 
liquids^^. Although we shall present results only for the 
2x2 plaquette degree of freedom, it can equally be ap- 
plied to the 4x4 superplaquette case. 



A. General formulation 

Let us return to the original Hamiltonian and ex- 
plicitly separate out the condensate mode in the opera- 
tors 7ia: 



ga + Pi, 



(22) 



The condensation will occur in a certain superposition 
of the singie-plaquette states. The real-valued multiplet 
ga plays the role of a condensate wavefunction (CWF)— . 
Here it is chosen to be spatially homogeneous, but in- 
homogeneous phases can also be included. The CWF is 
normalized to the condensate fraction: 



E9l = 



no- 



The non-condensate bosonic operators /3ia describe fluc- 
tuation corrections to the HF solution. If they are ne- 
glected, we naturally return to the results of Sec. HI. 
It is important to observe, however, that the HF ground 
state corresponds to the Bose condensation on each lat- 
tice site, not only in the fc = mode. 

The superfluid mean-field approximation amounts to 
enforcing the Schwinger boson constraint on average: 



(23) 



neglecting fluctuations in the condensate channel, and 
retaining only terms quadratic in f3 in the Hamiltonian 



H = Na 



/ino 



(24) 



where we abbreviated {ij)a = ((*.?))) a-nd matrix 

elements of H^^^^ are given by Eq. ([5]). 

The CWF Qa is determined by the Gross-Pitaevskii 
equation, similar to the HF equation pT|) : 



(25) 



which dcflnes the chemical potential fi and guarantees the 
disappearance of linear terms in /3 from the Hamiltonian 
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FIG. 14: Ground state energy for the self-consistent solution 
(main panel) and after the first iteration (inset). The crit- 
ical point Jl^ becomes a A-point. The absence of points in 
the main panel around is due to bad convergence in the 
simulation. 
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FIG. 15: Condensate fraction for the self-consistent solution 
(main panel) and after the first iteration (inset). Notice the 
shift of quantum phase transition points J2^'^. 



of Eq. ([M)) . It is clear that ga{na = 1) = i?^ and ^{uq = 
1) = So- In other words, Eq. ((25)) reproduces the results 
of Sec. Ill, if no is forced to be unity. Naturally, the first 
line in Eq. ([24]) coincides in this limit (up to the chemical 
potential) with the expression p2|) for Eq. Quadratic 
terms in Eq. (|24p represent fluctuation corrections to 
the HF results and constitute the focus of our analysis 
below. 

The next step is to transform the quadratic part {H2) 



0.45 




J2/J1 

FIG. 16: Self-consistently computed staggered magnetization, 
Afatag, for J2 < and columnar magnetization along the x- 
direction, Mcoi(a;), for J2 > J2^. 

of the Hamiltonian in Eq. (|24p into momentum space: 

k.a k.a 
+ 4.ga^5a./3L'/fcai}, (26) 

where 8^ = (cos + cos fcy, 2 cos cos ky) and fe is de- 
fined within the plaquette Brillouin zone (i.e., there are 
A^n k- states) . This Hamiltonian can be diagonalized by 
the Bogoliubov's transformation: 

a 

«-fc.=EH^«/3,, + J, (27) 

a 

to a new set of bosonic operators afcj/, which represent 
quasiparticle excitations and annihilate the new ground 
state: akv\^o) = 0. Of course, only positive quasiparticle 
energies, labeled by v, have physical meaning, however, 
in order to obtain closure relations for the wavefunction 
(""fca' '^ka) (which is, obviously, even in fe), we need to in- 
clude zero-energy eigenvectors as well^. This complete- 
ness relation has the form, valid for all wavevectors: 

''^i'^ka'^kb ^ "^'fca^'fch) =Sab', 

V 

The amplitudes uj^(fc) and w^(fe) are determined from 
Bogoliubov's equations: 

Vi{k)ul^ + Ul{k)vl^ = - ^Ak)vl^. (29) 
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where we have introduced symmetric matrices: 

a 

+ 25^(i7rnt)::>a,5a.; (30) 

(7 

a 

It follows from Eq. ([2^ that at each k the quasiparticle 
amplitudes obey the orthogonality conditions: 



E 

a 
a 



0. 



(31) 



For any value of A; Bogoliubov's equations always 
have at least two zero eigenvalues, which correspond to 
the zero-norm eigenvector u = —v = g. This means 
that our case differs fundamentally from the canonical 
superfluid Bose gas: instead of having a macroscopic 
number of particle in one particular energy state, we ob- 
tain a macroscopic number (equal to Nq) of condensation 
modes, each containing less than one boson. 

The quasiparticle energy equals 2uji,{k) and the GSE, 
condensate fraction and spin polarization are expressed 
in terms of and v]^^ as: 



^1 
N 8 



no 



k.u.a 



(32) 



In this expression fe-summations are extended over the 
plaquette Brillouin zone and (^-summations over positive 
eigenvalues of Eq. , as indicated by the primes. 



B. Results for the symmetric plaquette covering 

The condensate fraction no should, in principle, be de- 
termined self-consistently. The approximation is reason- 
able if riQ ^ 1- However, close to the phase transition 
points this is not true, since fluctuations are very large 
in their vicinity. But deeply in each phase the approxima- 
tion works reasonably well, because no turns out to be of 
the order of 0.9. Results of the numerical solution of Eqs. 
(PS)) . and ([5^ for the symmetric covering of the lat- 
tice with 2x2 plaquettes are shown in Figs. [TUfTTl The 
system size is 100 x 100 plaquettes and periodic boundary 
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FIG. 17: The lowest excitation branch uoiik) along the [10] 
and [11] directions for three values of J2/J1 chosen in different 
phases. 



conditions are assumed. Figures' main panels correspond 
to the self-consistent solution and their insets give results 
after the first iteration, which is equivalent to solving 
the time-dependent Gross-Pitaevskii equation^. Due to 
bad convergence close to the transition points (see, for 
instance. Fig. [Ti)) the values of and J^^ were deter- 
mined by extrapolation: w 0.33Ji and k, 0.65Ji. 
The large shift of compared to the HF value is due to 
fluctuations in the /3-channel, which renders this point to 
be a A-point, reduces the nominal value of the magneti- 
zation in the Neel phase down to Af ( J2 = 0) « 0.37 (Fig. 
I16p . and causes a great suppression of the condensate, as 
shown in Fig. [T5l 

However, the most interesting quantity to observe is 
the gap in the excitation spectrum. Due to the homo- 
geneity of the plaquette lattice, it occurs at fc = and is 
shown in Fig. [TTl 

Technically, one may show that its very existence re- 
flects the nature of the ground state in the paramagnetic 
phase. Indeed, introducing linear combinations of the 
amplitudes u and v: if = u + v and x = u — v, Bogoli- 
ubov's Eq. (|29p can be rewritten in the form: 



(C/^ + C/'^)(C/^-[/^)x 



In the non-magnetic phase the condensation occurs in the 
lowest plaquette state jllOO): ga = \/noSa,iioo and the 
chemical potential coincides with its energy: fj. = enoo. 
Moreover, the matrix (^int)h 1100 vanishes. Writing 
down the remaining matrices in (j30p at fc = 0, it is easy 
to see that there exists only one vector which is anni- 
hilated by ([/^ — U^) . Outside the intermediate region 
this simple situation is not valid and there exist three 
eigenvectors Xj which correspond to uj'^ = 0. One of 
them is the condensate mode and should be discarded. 
The other two give doubly degenerate Goldstone modes 
in the Neel and columnar phases. Here the self-consistent 
field, determined by ga, breakes the spin-rotational sym- 
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metry of the original Hamiltonian. However, since the 
CWF ga belongs to the M ~ subspace, the generator 
Sz remains an integral of motion. Thus, there should be 
two Goldstone modes associated with rotations around 
the X and y axes^. 

Our approximation correctly describes the excitation 
spectrum only at small fc. However, this is more than 



enough to observe that the collective modes are of the 
spin- wave type in the Neel and columnar phases, while in 
the paramagnetic phase the excitation band is parabolic. 
These conclusions are summarized in Fig. 1171 where we 
show the lowest branch uji{k) along two directions fc|| [10] 
and fc||[ll] for three values of J2/J1, chosen in different 
phases. 
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